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ABSTRACT 

We present a new automatic method to identify galaxy mergers using the morpholog- 
ical information contained in the residual images of galaxies after the subtraction of a 
smooth Sersic model. The removal of the bulk signal from the host galaxy light is done 
with the aim of detecting the much fainter and elusive minor mergers. The specific 
morphological parameters that are used in the merger diagnostic suggested here are 
the Residual Flux Fraction (RFF), and the asymmetry of the residuals (A(Res)). The 
new diagnostic has been calibrated and optimized so that the resulting merger sam- 
ple is very complete. However, the contamination by non-mergers is also high. If the 
same optimization method is adopted for combinations of other structural parameters 
such as the CAS system, the merger indicator we introduce yields merger samples of 
equal or higher statistical quality than the samples obtained through the use of other 
structural parameters. We investigate the ability of the method presented here to se- 
lect minor mergers by identifying a sample of visually classified mergers that would 
not have been picked up by the use of the CAS system, when using its usual limits. 
However, given the low prevalence of mergers among the general population of galax- 
ies and the optimization used here, we find that the merger diagnostic introduced in 
this work is best used as a negative merger test, i.e., it is very effective at selecting 
non-merging galaxies. In common with all the currently available automatic methods, 
the sample of merger candidates selected is heavily contaminated by non- mergers, and 
further steps are needed to produce a clean merger sample. This merger diagnostic has 
been developed using the HST/ACS F606W images of the A901/02 multiple cluster 
system(z = 0.165) obtained by the STAGES (Space Telescope A901/02 Galaxy Evo- 
lution Survey) team. In particular, we have focused on a mass and magnitude limited 
sample {log M/Mq > 9.0, i^vega, Total ^ 23.5mag) which includes 905 cluster galaxies 
and 655 field galaxies of all morphological types. 

Key words: methods: data analysis, galaxies: clusters: individual: A901/902; galax- 
ies: evolution; galaxies: interactions; galaxies: structure 



Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 



search in Astronomy, Inc. under NASA contract NAS 5-26555. 
These observations are associated with programme GO10395. 



2 Carlos Hoyos et al 



1 INTRODUCTION 

Mergers are the most extreme type of galaxy interaction, 
as the final product of a merger event can be totally dif- 
ferent from the original objects involved. Considerable ef- 
forts have been devoted towards the understanding of the 
physical processes that reg ulate galaxy me r gers starting 



close Dairs (IPatton et al. l200d I2OO2I; Lin et al.l 


2004; 


De ProDris et a 


[.1 12OO7I: Lin et al.1 12OO8I: lEllison et al.1 


2010; 


Robaina et al. 


2010|l, or even the correlation function 


(Beh et al.1 1200( 


J~Masiediet al.ll2006l). The moruholodcal 



pnysical processes tnat reg ulate galaxy me r gers startmg 
from the very early work of ISpitzer Baade (1 1951) or the 
seminal simulational works presented in Toomre Toonir3 
(|1972i ); IToomrd ((1973). These work made it clear that, 
even though the stars rarely collide with each other dur- 
ing a merger process, such episodes can have dramatic con- 
sequences for the gaseous component of the galaxies in- 
volved. Later works such as tho se of [B arnes & Hernquist 
(I1991L I1991 119961); iBarnesI (i200i) ; iBournaud et al.1 (l2005f); 



IWetzstein, Naab/fc Burkert (l2007l ) ; ISpringel HernquistI 



(I2OO5I) ; iBournaud et aD [20081) ; iHopkins et al.1 (|2009l ) 
IStewart et a l.' ('2009') ; 'Chihnga rian et al.1 (|2010') have helped 
address specific issues of merger processes such as the inter- 
nal structure of the remnants, the relevance of the orbital 
parameters or the impact of the gas fraction on the possible 
regeneration of galactic discs after a merger episode. 

However, the study of mergers is not only relevant 
because of the physics involved. The evolution of the 
massive early type galaxies that populate the red sequence 
cannot be explained using passive evolution models only, 
and mergers have been found to play a key role in their 
evolution. In particular, the evolution of the luminosity 
function and colours of galaxies since 2; ^ 1.0 observed 
in the COMBO-17 (Classifyin g Objects by Medium-Band 
Obse rvations in 17 Filters) (jWolf et al.1 l2003l ; iBell et al.1 
I2OO4I ) and the phase 2 of t he DEEP (Deep Extragalactic 
Evolutionary Probe) survey (|Faber et al.ll2007l ) suggest that 
the merger episodes have a huge impact on the evolution 
of early- type galaxies, increasing the stellar mass by a 
factor of two over the last 8Gyr. More re cently, the im- 
portan ce of mergers has been highlighted in iRobaina et al.l 
(|2010l ). who conclude that the evolution of massive, red 
galaxies depends strongly on their merging history. Further 
studies have helped ascertain the impact of mergers in 
specific aspects of the evolution of red gala xies, including 



mass es (Ivan Dokkum et aL 2010l), sizes (see Truiillo et al. 
20061. I2OO7I; iGiavaliscoTRavindranath. &; Daddil l2007l ; 
Buitrago etaP I2OO8I ; I van Dokkum et al.l l2008p among 



others) and velocity dispersions (ICenarro Truiilloi r2009) 

The identification of mergers in deep astronomical im- 
ages is thus a very important issue for galaxy evolution stud- 
ies, and the huge number of galaxies observed in modern 
surveys creates the need of reliable automated merger detec- 
tion mechanisms. A reliable merger identification technique 
is the key element in the calculation of the merger fraction. 
The merger fraction is defined as the fraction of galaxies 
with a recognizable ongoing merger episode that is found in 
any given (often mass-limited) sample. It is the first step 
towards the comoving merger rate, which is the number of 
merger events per Mpc^^Gyr""*^ . 

Several automatic identification techniques have been 
developed to single out mergers from non interacting galax- 
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20081; IJogee et al. 


|2009|), kinematical and spatial 



techniques are based on the fact that the objects involved in 
a merger episode will be gravitationally disturbed. The CAS 
(Concentration, Asymmetry, clumpineSs) system measures 
these specific aspects of the surface brightness distribution 
of galaxies in order to identify mergers. Objects with high 
asymmetries are usually taken to be mergers by this method. 
On the other hand, the G — M20 system measures whether 
the galaxies appear to be shredded or not, since both the 
Gini and the M20 numbers measure whether and how is the 
light concentrated in any given object. In this system, the 
more shredded galaxies are selected as mergers. The pairing 
techniques look for pairs of galaxies whose relative positions 
and velocities should be conducive to strong interactions in 
a relatively short timescale after observation. Each of these 
methodologies are sensitive to different time scales, mass ra- 
tios, orbital parameters, and gase o us con t ent of the galaxies 
involved. For instance, Conselice (20061); iLotz et al.1 (|2008l ) 
conclude that the CAS parameters are sensitive to roughly 
a timescale of 0.4 ^ 1.0 x 10^ years, while the time sen- 
sitivity of the pairing techniques depend on the projected 
separation between the galaxies. 

This work contributes to the morphological automated 
detection of mergers. Here, the morphological parameters of 
the residua l images afte r the subtraction of a smooth Sersic 
model (see ISersid Il963l , for a definition of this profile) are 
explored. This should, at least in principle, better reveal the 
impact of the gravitational interaction on the morphology 
of galaxies. This was done with the aspiration of detecting 
minor mergers. An isolated galaxy will, with time, adopt 
an approximately symmetric profile whereas an interacting 
galaxy will appear to be more asymmetric. The removal of 
an intrinsically symmetric profile such as the Sersic model, 
which could be regarded as the quiescent, underlying galaxy, 
will more clearly expose the asymmetric signature of the 
light from an interacting galaxy. Thus, the structure of the 
residuals is investigated with the aim of finding the combina- 
tion of structural parameters that produces merger samples 
of better statistical quality. This optimization step is done 
in an unambiguous way, by using an objective criterion to 
grade the performance of the diagnostics tried. The specific 
criterion used here encourages completeness at the expense 
of a fairly high contamination by non-mergers, and the re- 
sulting merger sample needs to be cleaned afterwards. 

In this paper, 31] presents the observational data used 
together with the galaxy samples selected to derive and then 
test the proposed method. Section JSJdescribes the data pro- 
cessing techniques employed, and the structural parameters 
used. Section SI Presents the objective method introduced 
here to determine what combination of structural parame- 
ters produces the merger sample of highest statistical qual- 
ity. The precise definition of the "statistical quality" of a 
sample that is selected from a parent population is given 
in ^4.11 Section presents a visual analysis of the objects 
selected as potential mergers by the method presented here, 
focusing on the contamination of the resulting sets of merger 
galaxies by non- mergers ( ^5.1|) and on the ability to detect 
minor mergers ( ^5.2|) . Finally, ^presents the conclusions of 
this work. 
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2 DATA 

The data used to illustrate this method is provided by the 
HST/ACS F606W observations that were obtained as a part 
of the Space Telescope A901/902 Galaxy Evolution Survey 
(STAGES)E1 (Gray et al. 2009). STAGES is a multiwave- 
length project that was designed to explore the impact of 
environment on galaxy evolution. Its main target is a mul- 
tiple cluster of galaxies located at z ^ 0.165 that harbors 
different environments with different densities. At the dis- 
tance of this multiple cluster, 1 arcsecond corresponds to 
2.83 kpc. This survey also includes X-ray XMM-Newton, 
UV GALEX, IR Spitzer, spectroscopic 2dF, Radio Giant 
Metrewave Telescope, and optical COMBO-17 observations. 

The HST/ACS STAGES observations form an 80-tile 
mosaic that covers almost 30 x 30 arcmin^ in the F606W 
filter, with an average exposure time of around 2 ks. The 
observations were reduced using an output pixel scale of 
0.03 arcsec, and a pixfrac of 0.8, in order to keep the PSF 
ellipticity as stable as possible for weak lensing studies. 
The PSF Full- Width-Half-Maximum (FWHM) is 3.12 pix- 
els. The point source completeness limit of these images is 
F606W,AB = 28.5mag. These observations are deep enough 
so as to probe most of the luminosity function of galaxies 
in clusters, as it is possible to recover reliable structural 
information for galaxies up to an absolute magnitude of 
Mf606W,ab = — IS.Omag. This limit is just two magnitudes 
brighter than the most luminous globular clusters, and is 
typical of dwarf elliptical systems. 

2.1 Sample Selection Morphologies 

This work makes use of a mass and magnitude selected sam- 
ple (9.0 ^ log M/M0, RveRa ^ 23.5mag) whic h is very simi- 
lar to the sample used in Malt bv et al T (|2010h . The galaxies 
are all in the 0.05 ^ 2;phot ^ 0.30 redshift range, with a 
relative dearth of sources in the lower end of this interval. 
The mass limit ensures that the sample is complete in stellar 
mass for both the blue clou d and the red sequence, as was 
shown in iBorch et al.1 (|2006l l. The magnitude limit guaran- 
tees reliable visual morphologies, since all the sources show 
extended images in the HST/ACS data. The sample includes 
1560 galaxies distributed among both the Field and the 
Cluster environments. The sample can also be divided into 
four different morphological classifications labeled as "E", 
"SO" , "Sp" , and "0th" . The "E" bin is made of course of el- 
liptical systems, the "SO" bin gathers the lenticular galaxies, 
the "Sp" bin comprises the spiral galaxies, and the "0th" 
bin includes a mixture of irregulars, compact, and highly dis- 
turbed sources that do not fit into any of the other galaxy 
classes. The "0th" bin includes 235 irregular galaxies, 41 
disturbed galaxies, and 13 compact sources. This bin is thus 
dominated by irregular systems. This sample was selected 
for the current work due to the following reasons. 

Redshift Uniformity: Most of the galaxies are found in 
the fairly narrow 0.16 ^ z ^ 0.30 range, and only 7% of the 
objects are in the 0.05 ^ z ^ 0.16 interval. The lookback 
time is thus approximately the same for most objects and 
no age related systematic effects are expected. The STAGES 

^ http:/ /www. nottingham.ac.uk/astronomy/stages 



HST/ACS observations will also probe similar rest frame 
wavelength intervals redwards of the 4000 A break, observ- 
ing the same stellar populations. The spatial resolution is 
approximately the same for most sources, too. There is no 
need to apply morphological k — corrections or distance cor- 
rections to the structural parameters. It is important to keep 
in mind, though, that the total /c— correction spread could be 
relatively large, ranging from 0.0 mags for the nearest blue 
objects to 1.0 mags for the farthest red galaxies. The {l-\-z)^ 
cosmological dimming could also introduce a bias of around 
0.5 mag between the high redshift objects and those located 
at the cluster's distance. However, the sample used in this 
work is a mass limited sample, and the faint magnitude cut 
applied is only made in order to ensure that all the relevant 
structural parameters can be measured reliably. This makes 
it possible to estimate the fraction of sources that could be 
affected by the k — correction and cosmological dimming bi- 
ases previously mentioned. This is done by adding 1.5 mag 
to the observed magnitudes and calculating the fraction of 
sources that would have been excluded from the sample in 
this situation. This simple estimate indicates that 90% of 
the galaxies in the current sample should be free from these 
systematics. 

Variety of Environments and Morphologies: This sample 
includes a good number of galaxies of several morphological 
types residing in very different environments. This makes 
it ideally suited to explore whether the performance of our 
merger indicators is sensitive to these variables or not. 

All the objects of this sample also have good morpho- 
logical information, which will be presented in detail in a 
forthcoming paper (Gray et al. in prep.). The morpholog- 
ical catalogue was created using visual classification from 
seven STAGES team members who were first trained on 
a consistent subsample of previously classified galaxies. All 
5090 galaxies with R < 23.5mag and 2;phot < 0.4 were clas- 
sified blind without knowledge of their cluster or field mem- 
bership. Galaxies were randomly assigned to classifiers who 
used a revis ed Hubble T-type scale and weighting scheme 
described in iLane et al.l (|2007[ ) . Each galaxy in the sample 
was classified by three people, while a subset of 786 bright 
galaxies previously studied in IWolf. Grav, Meisenheimerl 
((2005|) received classifications from all seven. 

In addition to the revised Hubble T-type, classifiers 
were also able to note further information about the galaxy 
structure within certain well-defined parameters. The dy- 
namical state involved an interpretation of the probable 
cause of any observed disruption, where the possibilities 
were a tidal interaction with a neighbour (I); tidal interac- 
tion suggesting a merger (M); tidal feature without obvious 
cause (T); and chaotic systems (C). The degree of disrup- 
tion was quantified by a disturbance parameter which was 
allowed values of 0, indicating little or no disturbance; [1,2] 
indicating moderate/strong asymmetry (e.g. an HII region); 
[3,4] showing a moderate or strong distortion (e.g. a tidal 
tail). 

The precise sample definition and its breakdown are 
now summarized. 

The Cluster Sample: This sample is defined by a redshift 
interval zphot = [0.17- A(i?), 0.17+A(i?)], where the photo- 
metric redshifts have been calculated from the COMBO-17 
data, and the half-width 
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A(R) = v^O.OlS^ + 0.00965252 x (1 + ioo-6x(RTot-2o.5)) 
is allowed to vary with apparent R band magnitude. This 
redshift interval varies with apparent brightness due to the 
higher precision of the COM BO- 17 photomet ric redshifts 
for the brighter sources, (see iGrav et al.1 l2009l . for a more 
detailed explanation.). The photo- 2; half- width distribution 
was normalized so the completeness was ^ 90% at any mag- 
nitude. For the faint end, there is also some sample contam- 
ination from field galaxies. This was calculate d using the 
counts of the smooth models seen in Fig. 14 of iGrav et all 
(|2QQ9l l. When estimating the contamination, the field distri- 
bution was assumed to be consistent with the average galaxy 
counts N{z, R) outside the cluster. The average contamina- 
tion is 30%. This subsample thus contains cluster galaxies 
with good photometry for which the HST/ACS observations 
show an extended source. This cluster sample contains 905 
galaxies. 

The Field Sample: This includes galaxies in the redshift 
intervals zphot = [0.05,0.14] and z = [0.22,0.30]. In this 
subsample, only HST/ACS extended sources with good pho- 
tometry are included. This field sample contains 655 galax- 
ies. 

This leaves a final sample of 655 Field galaxies (100 E, 
60 SO, 318 Sp, 177 0th) and 905 Cluster galaxies (192 E, 216 
SO, 383 Sp, 114 0th). As before, the "0th" galaxies represent 
a mixture of irregulars, compact, and highly disturbed sys- 
tems. This classification is dominated by irregular galaxies. 
The STAGES morphological catalogue has good morpho- 
logical information for all objects in this sample. The source 
detection runs presented in ^3.11 and the structural parame- 
ter analysis presented in ^3.2l vielded structural information 
for a total of 1537 sources. Thus, there is structural infor- 
mation available for 98.5% of the full sample. 



3 DATA PROCESSING 

This section presents the sextractorEI 
(|Bertin Arnout3 ri996) runs and subsequent mor- 
phological analysis carried out here. The main ingredients 
this work needs are an error or weight image and a mask 
image. The err or ima ge is needed by both SEXTRACTOR 
and GALFIT (Peng et al.t 12002 ) in order to modulate the 
importance of each pixel in the photometry and morpho- 
logical properties. The mask image, which is created by 
SEXTRACTOR, is later needed by GALFIT so that it 
fits galaxies in a way that is designed to identify mergers. 
The basic idea is to produce a simple Sersic model that 
is as focused on the main body of the target galaxies as 
possible. The putative mergers should not therefore have 
any effect on these fits and will be either simultaneously fit 
or masked out. 

This Sersic model should be an appropriate description 
of the surface brightness distribution of elliptical galaxies. 
However, it has to be noted that in the case of lenticu- 
lar and spiral galaxies, bulge+disc decompositions (B+D) 
would produce better descriptions. However, there are a 
number of reasons that make simple Sersic fits preferable. 
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• A simple Sersic profile will remove a fraction of the 
smooth, symmetric signal from any given galaxy, regard- 
less of its morphology. This statement remains true for disc 
galaxies. 

• The goal of the current work is to develop an automatic 
and fork free method. The use of B+D decompositions to- 
gether with simple Sersic fits would require an additional 
step in order to decide what functional form describes best 
each galaxy. There are objective ways to do this , such as 
the Bayesian method presented in lden Brok et aP (|201lh or 
the Akaike information criterion, but this is well beyond the 
scope of the current work. 

• B+D decompositions need to be re-examined for their 
internal consistency, as their final results might depend on 
the initial point in parameter space chosen for the galfit 
minimization. 

Thus, it is not possible to us e the SEXTRACT OR and 
structural catalogue presented in lCrav et al.l (|2009l ) because 
the source detection and fits shown there were made with the 
goal of creating a robust catalogue of structural parameters 
derived from single Sersic fits. Instead, the SEXTRACTOR 
configuration used in the current work (see ^3.ip is fairly 
sensitive to the detection of faint, small features, as it is 
set to detect the faintest sources that can be said to be de- 
tectable in the images used. More importantly, the deblend- 
ing parameters are selected so that SEXTRACTOR will 
tell apart sources with up to a three magnitude brightness 
difference, incorporating much of the 1 : 3 to 1:10 mass ratio 
range defining minor mergers, assuming that their M/L ra- 
tios are similar. These SEXTRACTOR parameters are thus 
not only intended to separate the objects involved in a major 
merger, with luminosity ratios between 1 : 1 and 1:3, but are 
also geared towards the deblending of the objects involved in 
a minor merger episode, with luminosity ratios between 1 : 3 
and 1 : 10. Also, the very sensitive detection threshold em- 
ployed makes it more likely that the less luminous segment 
of any given deblended object will be engulfed by the Kron 
aperture of the larger segment. This is key in the current 
analysis, since it is hypothesized here that such less lumi- 
nous segments could be the less luminous galaxies in merger 
episodes. However, these smaller segments could just be HII 
regions or simply objects along the same line of sight, which 
will naturally lead to a contamination of the merger sam- 
ple. This will be discussed in ^5.11 The larger segments in 
deblended sources woul d then be th e more luminous galax- 
ies. The Kron aperture (|Kronlll98Ql ) is defined to contain a 
specific fraction of the light of a galaxy. Such fraction de- 
pends on the intrinsic profile of galaxies, but it ranges from 
90% for the steepest profiles to 95% for exponential disks. 
Both percentages depend on the observed surface brightness, 
though, in the sense that SEX TRACTOR misses more flux 
from the dimmest objects. See iHammer^ t al. (201^) for a 
description of this systematic error. 

The larger and more luminous object is then fit by a 
smooth Sersic model which is created using GALFIT ac- 
cording to the rules explained in ^3.21 This model is then 
subtracted from the original image in order to estimate how 
the image would look without the more luminous galaxy 
of a merger event. This residual image still contains most 
of the signal from the less luminous sources that were found 
within the observed Kron aperture of the larger object. This 
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ensures that these smaller sources will have a great impact 
on the structural properties of the residual images. This is 
indeed one interesting property of these residual images that 
will be exploited. 

The structural parameters of both the original image 
and of the residual image are then calculated within the 
aforementioned Kron aperture. The specific morphological 
parameters used are shown in ^3.31 These are mostly based 
in previous works. The structural parameters of the residual 
image within the Kron aperture are then expected to be 
very sensitive to the smaller and less luminous member of 
the merger. 



3.1 SEXTRACTOR Runs 

SEXTRACTOR was run twice on each of the 80 tiles that 
compose the STAGES HST/ACS F606W mosaic. The first 
pass was performed for the sole purpose of obtaining an em- 
pirical map of the background variance and of the filtered 
poissonian signal. Such images can be produced by SEX- 
TRACTOR as output images. These images were then com- 
bined to produce an appropriate weight image to be used in 
the second pass. The first SEXTRACTOR run is presented 
in ^3.1.1[ and it can be safely skipped by users providing 
their own error or weight images. The second run is pre- 
sented in 33X21 which presents the SEXTRACTOR con- 
figurations used to produce the final source lists and mask 
images. 



3.1.1 Preliminary SEXTRACTOR run 

The most relevant SEXTRACTOR parameters used in the 
first run are summarized in table [T] Table [1] also presents the 
main results obtained from this first SEXTRACTOR run, 
averaged over the 80 tiles. The average background sigma 
is an estimate of the noise that is observed in the image 
areas in which there are no galaxies. This is set by the read- 
out noise of the Analog-to-Digital converters of the ACS 
camera, the poissonian noise of the sky background and the 
use of the MULTIDRIZZLE technique. It was estimated as 
the average of the median values of the background RMS 
images created by SEXTRACTOR in this first pass. The 
average effective gain (Geff) roughly measures the growth 
in the photometric errors caused by the intrinsic poissonian 
nature of photon counting measurements. It was measured 
using the background subtracted filtered frame and a very 
rough empirical estimate of the per-pixel RMS image ob- 
tained b y means of a loose ad aptation of the method pre- 
sented in iGrazian et al.l (|2006l ) , which shows precise formu- 
lae to calculate the RMS of an image in the case of noise 
correlation. The effective gain image is then. 

where I is the background subtracted filtered image, and 
RMS is the empirically derived uncertainty. 

The (Ceff) number is yielded by the value of this im- 
age in the brighter areas of the images dominated by the 
poissonian noise. 

The measured value of (Ceff), reported in table [1] agrees 
very well with its expected value of around 1445, which can 
be theoretically estimated as: 



(2) 



where GeS is the effective gain, g is the original detector's 
gain (2, for the STAGES observations), T is the total ex- 
posure time, and the fraction is the ratio between the ef- 
fective areas of the pixels before and after the operation of 
the MULTIDRIZZLE technique. The latter value of 1445.0 
was adopted for use in the second run of SEXTRACTOR. 
The measured value of Ceff was then merely used as a san- 
ity check to ensure that the errors in the input image indeed 
behave as expected. 

The first SEXTRACTOR run thus produces an esti- 
mate of the background aekg and a background subtracted 
filtered frame containing the poissonian signal S. These were 
combined according to the usual CCD error equation in or- 
der to obtain a weight image for use in the second SEX- 
TRACTOR run: 



Weight = 



'Bkg 



+ S/Ge 



(3) 



where Weight is the final weight image, as kg is the back- 
ground RMS image, S is the background subtracted filtered 
image created by SEXTRACTOR, and Ges = 1445.0. The 
final weight image was later processed using the WEIGHT- 
WATCHER^ code to ensure that problematic pixels with 
either zero exposure time (typical of the image edges) or 
with saturated signal were assigned zero weight. Less than 
2% of the pixels had to be discarded in this final procedure. 
WEIGHTWATCHERS also produces a flag image that was 
used in the second pass of SEXTRACTOR. This final 
weight frame is then used in the second run of SEXTRAC- 
TOR to supress the detection of objects in low weight pixels 
and give appropriate pixel weights for the photometry. 



3.1.2 Second SEXTRACTOR run 

This second pass was used to obtain the final segmenta- 
tion image and a SEXTRACTOR catalogue, which are key 
to the GALFIT analysis of the images. The segmentation 
images separate object pixels from background pixels, and 
provide basic photometric information to be used as initial 
conditions for GALFIT. It is highlighted here that the in- 
put photometric catalogue of target sources is not produced 
by this SEXTRACTOR run. The list of objects to fit and 
study is defined in and SEXTRACTOR is run here 

with the sole purpose of obtaining basic photometric infor- 
mation about this pre-defined sample. The new configura- 
tion used is presented in Table [2] 

As it can be seen from Tabled the new configuration is 
fairly aggressive . It is more sensitive than the "hot" config- 
uration used in iGrav et al.1 (l2009l ). The minimum nominal 
integrated S/N of the detections is 8.3. This is 1.7 times 
higher than the usual S/N limit of 5.0 which is usually ac- 
cepted for a point source detection. The new SEXTRAC- 
TOR parameters were chosen this way because, to zero- 
th order, the effect of noise correlation is an artificial in- 
crease of the a image with respect to the inverse square 
root of the weight image. Thus, a stronger signal is needed 
in order to spawn a genuine detection. Given the number 



See http : //www . astromatic . net/sof tware/weightwat cher . 
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Table 1. SEXTRACTOR parameters used for the first pass of this code on the 
HST/AGS F606W images. This first run was used to obtain an empirical estimate of 
the final error image. The main results obtained from this first pass are also presented 
in this table. 



of images that were multidrizzled together and the target 
resolution of the HST/ACS images used, experience shows 
that the nomin al S/N has to be multiplied by around 0.7, 
as is shown in iHammer et al.l (|201 0) and Case rtano et al.l 
(2000 ). This is the reason behind the additional factor of 
1.7. The second SEXTRACTOR configuration makes use 
of a gaussian filter with a FWHM of 3 pixels, and a min- 
imum detection area of 12 pixels. This effectively removes 
spurious noise peaks, which are generally much smaller than 
the instrumental FWHM. As for the deblending parameters, 
the ones shown in Table [2] would, in principle, make SEX- 
TRACTOR separate objects with a factor of 100 difference 
in flux. However, the value of the DEBLEND_NTHRESH param- 
eters also affects the deblending process. With the choice 
of 32 deblending thresholds, the child detections tend to 
bear a larger fraction of the total flux than just 1% of the 
parent source. This was done this way with the idea of sep- 
arating minor mergers in a very late stage, with a mass 
ratio close to 1 : 10. The SEXTRACTOR parameters used 
in the second pass also ensure that the source catalogue 
will include detections all the way up the real, almost point 
source detection limit of the images. At the same time, the 
minimum nominal S/N ratio required to foster a detection 
and the DETECT_MINAREA value adopted ensure that the num- 
ber of spurious sources caused by the multidrizzle algorithm 
and noise peaks will be kept at a minimum. This configura- 
tion, however, is open to the inclusion of spurious extended 
sources caused by statistical fluctuations. Also, this configu- 
ration makes SEXTRACTOR include pixels with poor sig- 
nal within the isophotal area of the detections which guar- 
antees that the sky portion of the segmentation image is 
free from most of the flux originating from the uncovered 
sources. This also has the important effect of making the 
Kron apertures as large as they could possibly be given the 
data used. This maximizes the probability of a small galaxy 
in a minor merger event falling within the Kron aperture 
of the more luminous object. This is particularly important 
since the structural parameters measurements shown in ^3.31 
are performed over this Kron aperture. On average, SEX- 
TRACTOR found 34600 objects in each of the frames. This 
number i s to be compared w ith the total number of sources 
found bv lGrav et ahf (|2009h . which is 75805 in all 80 tiles. 
This is explained by the fact t h at the SEXTRACTOR con- 
figuration used in iGrav et all (|2QQ9h was optimized to find 
and fit i?ap ^ 24.0mag counterparts from a previous cata- 
logue obtained from the R — band COMBO- 17 data, while 
the goals of the SEXTRACTOR catalog used in this work 
encourage the detection of faint sources near the brighter 



ones. This does not mean that all the sources found by this 
second SEXTRACTOR run are legitimate, bona fide detec- 
tions. In fact, specific simulations indicate that the majority 
of objects with MAG_ISO dimmer than 27.0 mag are spurious 
objects caused by statistical fluctuations, while the majority 
of the objects with MAG_ISO brighter than 27.0 mag are real 
sources. This is not a problem for the target sources studied 
in this work, whose MAG_ISO are all brighter than 24.0 mag. 
The fitting scheme used here, which is presented in ^3.21 just 
discards these faintest detections and therefore they will not 
have any impact on the actual fits produced for the much 
brighter targets of interest. At the same time, this SEX- 
TRACTOR configuration will detect and more importantly, 
isolate sources in the 24.0mag < FQOQW{AB) < 27.0mag 
magnitude interval, which are much more likely to be real 
objects. The comparison between the num ber of detections 
obtained here and in the Grav et al] (|200^ 1 work merely re- 
flects that the configuration used in the current work is much 
more sensitive to the smallest features, which increases the 
chance of including spurious detections in the SEXTRAC- 
TOR catalog. 

3.2 Galaxy Fitting: The GalP-Hyt Wrap Scripts 

The main contribution of the current work is that it begins 
to explore whether the morphological information contained 
in the residual images of galaxies can be used to assess if a 
galaxy is involved in a merger episode. Such residual im- 
ages are created by subtracting smooth models of the target 
galaxies from the original images while leaving most of the 
signal from other, possibly interacting objects in these resid- 
ual images. There are a number of codes capable of produc- 
ing models of galaxies in astronomical images by perform- 
ing two-dimensional model fits to their surface brightness 
distributions. The most commonly used ones are GIM2D 
( ISimardllTooil ) and GALFIT. Even though these codes are 
different in their specific details, their basic principles are 
very similar. They both try to minimize a possibly weighted 
value that depends on the structural parameters of the 
object being fit. This residual sum of squares is formed from 
the difference between the observational data and a trial 
function that is created according to a user-supplied set of 
rules. Chief among these rules are an error image, a mask 
image, and a PSF image. The error image regulates the rela- 
tive weight that the different pixels should be given. Thanks 
to this image, it is possible to prevent saturated pixels from 
having any weight in the figure- of- merit that the codes are 
set out to minimize. This image can also be used to reduce 
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Table 2. SEXTRACTOR parameters used for the second pass of this code on the 
HST/ACS F606W images. This second pass makes use of the empirical error image 
created using the output images of the first SEXTRACTOR run. The flag image 
created by WEIGHTWATCHERS was also used to exclude pixels with saturated 
signal or zero effective exposure time. 



the impact of areas with lower or no exposure time in mo- 
saic images created using the MULTIDRIZZLE technique. 
The mask image can be used to modify how the different 
fitting codes treat the different areas of the input image. 
For instance, it can help the fitting codes to tell what pixels 
belong to the target object being analysed, and what pixels 
belong to other objects and should therefore be discarded. 
The PSF image is also a key ingredient in the structural 
analysis of the surface brightness distribution of galaxies. It 
is also possible to constrain the parameter space region that 
the codes are allowed to explore. 

The present work uses the GALFIT code. The setup 
with which GALFIT runs is created by a python code called 
GALP-HYtQ written by CH. The GALFIT setup used is 
designed to create and subtract a smooth model of the pri- 
mary or target sources, while leaving the signal of nearby ob- 
jects largely intact in the residual images. In this context, the 
target sources are the galaxies that are being studied, and 
close sources are the objects found within the Kron aperture 
of the primary galaxy. Given the SEXTRACTOR param- 
eters used, these close sources could be the less luminous 
galaxies in a minor merger event, or one of the protagonists 
in a major merger episode. Thus, GALFIT is configured to 
produce a residual image in which only the model for the 
primary source has been removed. This is done in order to 
guarantee that the Kron aperture in the residual images will 
contain most of the information from the merged sources, 
while being free from the effect of the target galaxy. The 
GALP-HYT code is now briefly described. 

For each science image or tile it is given, it takes the 
following information as input: 

• The science image itself. These are the same im- 
a^es used in the structural parameter catalog presented in 
iGrav et al.l (|2009h . 



^ Pronounced Galp-Hit. 



• The weight image that SEXTRACTOR used during its 
second pass. This is converted into a a— image by taking its 
inverse square root. In this process, pixels with zero weight 
are given a very large value of a. 

• The SEXTRACTOR list of sources and segmentation 
image created during the second pass of the detection soft- 
ware. 

• A list of primary targets. This list of targets is given in 
a separate text file, one object per line. In this file, targets 
are identified by their SEXTRACTOR number IDs. 

• A PSF imagelet. Th is was again taken from the work 
done in lGrav et al.l (|2009[ ). This is a PSF imagelet that was 
built from many different non saturated stars found across 
a number of tiles and hence its S/N ratio is very high. 

For each object in the list of targets the SEXTRAC- 
TOR catalogue is examined and objects are classified with 
respect to the primary object according to the following set 
of rules: 

• The target object itselL It is fit using a single Sersic 
model with a free floating diskiness-boxiness "CO" param- 
eter. This extra freedom allows the model for the target 
object to take into account a larger fraction of the symmet- 
ric, undisturbed signal of the primary source. At the same 
time, this extra parameter does not complicate the interpre- 
tation of the fits since it is unlikely to introduce important 
degeneracies. Also, the evolution of "CO" as a function of 
the number of minor merger events is a prediction of the 
iBournaud et al.l (|2005l . l2008f ) models. 

• Objects up to two magnitudes fainter than the target 
source whose centres lie within the Kron aperture of the tar- 
get source are tagged as "A". In this context, the relevant 
magnitudes are the MAG_ISO magnitudes. The two magnitude 
difference was used since this is approximately the expected 
magnitude difference between the galaxies involved in a mi- 
nor merger with a mass ratio of 10:1, assuming that the 
M/L ratios are similar. These sources are fit with a single 
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Sersic profile with elliptical isophotes. The centres of the "A" 
components are fixed to the values found by S EXTRAC- 
TOR. This constraint prevents degenerate fits, as the "A" 
components will not blend and shift towards the primary 
object. 

• Objects more than two magnitudes fainter than the tar- 
get source whose centres lie within the Kron aperture of 
the target source are tagged as "B" . Again, the magnitudes 
used in this criterion are the MAG_ISO magnitudes. These 
sources are fit as simple exponential models with elliptical 
isophotes. This ansatz is less flexible than the one used for 
the "A" sources, but this is justified by the fact that there 
is less information available for these dimmer sources. The 
use of exponential profiles also helps to increase the execu- 
tion speed. Again, the centres of the "B" components are 
directly taken from the values found by SEXTRACTOR in 
order to prevent degenerate fits. 

• Objects up to three magnitudes fainter than the target 
source outside the Kron aperture of the target source are 
tagged as "E" (for external). As before, MAG_ISO measure- 
ments were used. These objects are fit using simple exponen- 
tial disks. Objects more than three magnitudes fainter than 
the source of interest are ignored. The centres, ellipticities 
and position angles of the "E" components are again taken 
from the SEXTRACTOR catalogue. 

The segmentation image produced by SEXTRACTOR 
is also modified so that the pixels belonging to the target 
source and to the "A" objects are given a value of 0. This 
ensures that these pixels are taken into account by GALFIT 
when calculating its figure-of- merit. Pixels belonging to "B" 
and "E" objects are not nulled, and thus the pixels contained 
within their ISOAREAS have no weight in the fit. Only the ex- 
tended tails of the "B" and "E" objects that go beyond their 
ISOAREAS are fit by the exponential profiles. It is important 
to note that there will be many objects three or more mag- 
nitudes fainter than the target source that will not fit any 
of the three categories above. These faint objects are not fit 
in any way and have no effect on the fits. This is important 
because the integrated magnitudes of the sources of interest 
are all brighter than F606W{AB) = 24.0mag, and therefore 
the large number of spurious sources with magnitudes fainter 
than F606W{AB) = 27.0mag wih have no effect on the fits. 
The sky is left as a free floating parameter. This is justified 
by a number of reasons. First of all, the SEXTRACTOR 
configuration used guarantees that the isophotal apertures 
of the detected objects reach to fairly faint surface bright- 
nesses, and thus the sky pixels in the segmentation image 
contain very little residual signal from the detected objects. 
Also, the masking process implies that the majority of the 
flux from the "B" and "E" objects is not taken into account 
by the fits and only their tails are fit using exponential pro- 
files. This naturally means that whatever their contribution 
might be to the average sky level affecting the target object, 
it will be approximately corrected for by the exponential fits. 
Finally, the "A" objects are fully fit. For these reasons, it is 
appropriate to leave the sky as a free fioating parameter. 

The size of the fitting box is defined in two steps. In the 
first step, the maximum between 150 pixels and the circu- 
larized Kron diameter is used. In the second and final step, 
this area is expanded so that it encompasses the centres of 



all "E" objects whose Kron apertures intersected the first 
box. 

The total number of additional sources that have to be 
fit in conjunction with a primary target is around 20. It is 
highlighted here that each SEXTRACTOR detection is fit 
by a simple, solid profile. In each case, the target object 
is fit by a single Sersic model, and the majority of the re- 
maining components corresponds to "E" objects and do not 
merge with the main body of the target galax y. This number 
of ext ra components is not unusual, see e.g. IHaussler et al.l 
(|2007l ) . The main target galaxy is fit using a single Sersic 
model. In most cases, the majority of the remaining compo- 
nents are external (E) objects whose main purpose is to help 
galfit compute a good sky value. The A and B components 
are the putative minor mergers, and their role is to allow 
the Sersic profile to provide a good fit to the main body of 
the target galaxy. 

Although this is very expensive in terms of CPU time, 
the GALP-HYT code was written to be run in High Per- 
formance Computers with thousands of CPUs, as it features 
a semi- intelligent built-in system of organizing its internal 
data fiow. It was run in the HPC computer of the Univer- 
sity of Nottingharr[3- Initial values and constraints are taken 
from the source catalog created by SEXTRACTOR. 

3.3 Additional Structural Parameters. The CAS 
Indices and the Gini-M2o System 

One of the common ways to tackle the automated detection 
of mergers is based on the signatures that mer ger events 
leave o n the morphology of ga laxies. As shown in IConselicd 
(|2003h : IPapovich et al. (|2005t ), the structures of galaxies 
contain important information about their past star forma- 
tion modes, and they can also shed light on their interaction 
history. These ideas are used to identify mergers according 
to their morphological properties. In the current work, the 
structural properties of the residual images after the sub- 
traction of the Sersic model described in 13.21 are explored 
and used as merger diagnostics. 

There are two systems which are currently in use 
for the morphological identification of mergers. The 
first one is the use of the Concentration, Asymmetry, 
Clumpi neSs (CAS) num bers, which was introduced in 
iBershadv. Jangren. &; Consehcd (|2000h : IConsehcd (120031 ). 
This system has been used in many works aimed at 
the study of the merge r fraction in many different con - 
texts, see for inst ance Conselice. Yang, &; Bluckl (|2009[ ), 
Lopez-Sanj uan et al. (2009), and Jogee et al.l (|2009|). 

This system makes use of three different indices: 

• C: The concentration index "C" measures to what extent 
the light in the galaxy is concentrated towards its centre. It 
is defined as: 

C = 5.0 X log ( — ) (4) 

\r2o/ 

where rgo is the circular radius containing 80% of the total 
light from the galaxy, and r2o is the radius of the circular 
aperture that encloses 20% of the total light of the target 
galaxy. The concentration index "C" takes values between 

^ See http : //www . nottinghcim . ac . uk/hpc/ 
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2.2 and 5.0. If it is calculated for simple Sersic models it 
depends mainly on the Sersic index and then on the ellip- 
ticity. 

• A: The asymmetry "A" measures to what extent any 
given image changes under a 180 degree rotation around the 
point that minimizes the asymmetry of that image. It is 
defined as: 



Bl 



(5) 



where Uj represents the original image, and llj^ is a 180 de- 
gree rotated version of the original image. In the same man- 
ner, Bij is a patch of background, and B}j^ is a 180 degree 
rotated version of this patch of background. This contribu- 
tion from the background is minimized independently in the 
same manner. This asymmetry measurement is defined even 
for images whose average value is 0.0, as it is normalized to 
the sum of the absolute values of the fluxes from each pixel. 
The rotation centre is optimized so that the value of the 
first term in the subtraction is a minimum. Here, the rota- 
tion centre is allowed to lie at most 9 pixels away from the 
SEXTRACTOR defined centre. However, when calculating 
the assymetry of the residual images, the rotation center 
is only allowed to move 4 pixels. The second term in the 
subtraction also undergoes this optimization process and it 
removes the contribution to the asymmetry from the back- 
ground. The asymmetry index of real images of galaxies can 
take values between 0.0 and approximately 0.8. Most objects 
have asymmetries lower than 0.2, though. The asymmetry 
index of the residual images of galaxies after the subtraction 
of a single Sersic model ranges between 0.4 and 1.6. 

• S: The dumpiness "S" quantifies the fraction of light in 
a galaxy that is contained in clumpy distributions. Large 
values of S imply that the light if the galaxy is accumulated 
in few, distinct structures. Low values of S indicate that the 
light distribution is smooth. It is defined as: 



S= lOx 



T- ■ 



where lij again represents the original image, and I^j repre- 
sents a blurred version of it, which is produced by convolving 
the original image with a two-dimensional circular gaussian 
kernel with a typical dispersion of a. It is usually correlated 
with the size of the target galaxy. The residual image after 
this subtraction only includes signal that is included in high 
frequency features of the galaxy. Also, the convolution pro- 
cedure is applied to a blank patch of sky in the image. This 
ensures that the contribution from the background noise is 
discounted from the final value of S. The "S" parameter can 
take values between -0.5 and 1.5, although it depends on the 
size of the convolution kernel used. 

The second system is based on the use of the Gini 
an d M20 parameters, which were origin ally introduced 
by [Abraham, van den Bergh, Naid (|2003h and lLotz et alj 
The Gini coefficient G measures the light concen- 
tration, like the C parameter, but it is insensitive to any 
particular centre. It is calculated according to the following 
formula: 



G 



1 



I/I X n X (n - 1) ^ 



(2xi-n-l)x|/,| 



(7) 



where n is the number of pixels, fi is the flux observed in 
the i^^ resolution element, and the sum is made in ascending 
order of fluxes, so that fi-i ^ fi ^ fi+i- This G index 
tells whether the light is evenly distributed among the dif- 
ferent resolution elements of an image. The G index has a 
value of 0.0 for flat light distributions, and it has a value of 
1.0 for light distributions in which all the light is contained 
in a single pixel. In practical terms, the Gini coefficient of 
real galaxy images lies in the [0.35,0.85] interval. 

The M20 parameter is based on the second-order mo- 
ment of the light distribution Mxot- It is defined as: 

^20=log I ^^l-i^A __ I (8) 



where A is the aperture within which this number is ob- 
tained, fi is the flux of the i^^ resolution element, Lxot 
is the total apparent luminosity contained in the aperture 
used, Xc and i/c are the coordinates of the bary centre of the 
light distribution for which the index is being calculated, and 
Mxot = EiGA^^ ^ ((^-^c)^ + (2/-2/c)^). In this definition, 
and in contrast with the definition of G, fi-i ^ fi ^ fi+i- 

M20 measures how far from the galaxy centre it is pos- 
sible to find the brightest features of the surface brightness 
distribution of the light. The M20 number can go from -3.0 
for very concentrated objects to -0.4 for objects with shred- 
ded light distributions. 

In addition to the aforementioned indices, the residual 
flux fraction {RFF, see Hovos et al. 2011) is used. It is here 
defined as: 



RFF = 



EjGALFIT 



(9) 



where A is the particular aperture used to calculate this 
index, and /^alfit ^^le model created by GALFIT. 

The RFF as defined here measures the fraction of the 
signal contained in the residual image that can not be ex- 
plained by fluctuations of the background. The 0.8 factor 
included in the definition of the RFF ensures that the ex- 
pectation value of the RFF of a purely gaussian noise error 
image of constant variance (as opposed to a spatially varying 
variance) is 0.0. This fact arises from the following integral: 



0.^ 



27r / 

J —c 



dx. 



(10) 



which calculates the expectation value of the absolute value 
of a gaussian random variable. 

In the current work, the structural indices presented 
above were calculated for three different images. These are 
the imagelet containing the target object cropped by GAL- 
FIT, the simple Sersic model created by GALFIT and the 
residual image obtained by subtracting the second image 
from the first one. It is stressed that this latter image is 
only stripped from the signal of the main body of the tar- 
get galaxy. The indices were also calculated for an artificial 
image simulating the background noise that affects the pri- 
mary galaxy, which will be described later on. This image 
is needed to calculate the background terms in the A and 
S indices, but it is also used to estimate the errors in the 
derived morphological parameters. 
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In all cases, the aperture used was the Kron aperture 
calculated by SEXTRACTOR. This aperture was chosen 
for the following reasons: 

• This aperture is designed to trace an elliptical region in 
the image in which the contribution in flux from the source 
of interest to which it is associated is either dominant or 
noticeable. Outside the Kron aperture there is still flux from 
the primary source, but it is heavily affected by noise. It is 
therefore pointless to go much further since morphological 
perturbations of the target galaxy at these levels will be 
impossible to measure reliably. 

• The Kron aperture misses a definite fraction of the light 
from most common profiles. This fraction depends on the 
intrinsic shape or profile of the target source and on its ef- 
fective surface brightness as well. Thus, the Kron aperture is 
di S/N matched aperture that grabs a more or less constant 
fraction of the total flux for objects of similar S/N, with a 
slight dependence on its profile. 

• This aperture has a radius that is typically 60% larger 
than the Petrosian radius defined in Petrosian (1976 ). This 
is a little bit smaller than the typical aperture of choice in 
which the structural parameters defined above are calcu- 
lated in other studies, which is twice the Petrosian radius. 
The choice of the Kron aperture will therefore provide higher 
S/N measurements of the morphological parameters, losing 
only a minimal amount of information about the outer struc- 
ture of the studied objects. 

The whole Kron aperture is used to calculate the CAS, 
G, and M20 numbers for the real image, the model and 
the background noise frame. However, when calculating the 
structural parameters for the residual image, a small area 3 
pixels in diameter is removed from the centre of the Kron 
ellipse so that the indices are not biased by the fit uncer- 
tainties, PSF mismatches, and resampling problems in these 
complicated regions. This latter aperture in which the centre 
is excluded is also used to calculate the RFF and S struc- 
tural parameters. Also, the rotation axis for the calculation 
of the asymmetry parameter of the residual images is not al- 
lowed to drift from the optimal centre found for the original 
image by more than 4 pixels. In all cases, the background 
was subtracted from the images using the sky value yielded 
by GALFIT. This was done with the purpose of minimizing 
the impact of the background terms in the A and S numbers. 

For the calculation of the S parameter, a gaussian kernel 
with cr = 0.2 X Rk was used. In this expression Rr is the 
radius of the Kron aperture calculated by SEXTRACTOR. 

The background noise image that is needed in order 
to calculate the background terms of the A and S numbers 
was created from the RMS image that GALFIT used. It was 
created on an object by object basis, so that each object has 
an individualized background noise image. The first step is 
to change the data number values of the pixels in the RMS 
image that, according to the segmentation image created by 
SEXTRACTOR, have been flagged as objectETpixels. Their 
new value is then set to the median value of the remaining 
sky pixels in the RMS image. This modified image is then 
multiplied by a white noise image with cr = 1.0. This final 
step creates an image that is a good representation of the 
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Table 3. Typical absolute errors in the structural parameters 
used in this work. The left columns present the typical errors for 
the structural parameters of the original galaxy image, while the 
columns to the right show the uncertainities for the structural 
parameters of the residuals. Errors above the horizontal line are 
the random errors, while the entry below the horizontal line is 
the systematic uncertainties for the A(Obj) index. 

underlying noise that affects the measurements. It would be 
ideal if this image included correlated noise, but since the 
target galaxies are all larger than the error correlation length 
this image was deemed sufficient. 

The noise image was also used to calculate the errors 
in the structural parameters. This was done by using this 
frame to recreate ten realizations of the original image. The 
structural parameters were then recalulated for this set of 
realizations and an error estimate is obtained by the very 
simple prescription of removing the smallest one and the 
largest one. According to Chebyshev's inequality this pro- 
duces an interval whose upper limit is 2.3 3.0 x a. The 
actual value depends on the assumed underlying distribu- 
tion. Here, the error distribution is assumed to be gaussian 
for all structural parameters, and robust estimates of the 
1 — a uncertainty are obtained by doubling the Chebyshev 
error estimate. 

The first term of the A index was also calculated for the 
GALFIT model. This was used to estimate the systematic 
uncertainties in the asymmetries of the real galaxies, since 
they should be zero in the model image. Deviations from zero 
thus reflect inaccuracies in the minimizing algorithm. These 
deviations amount to less than 0.05 in almost all cases. 

Table [3] presents the typical errors in all the derived pa- 
rameters. Although each object should have its own error for 
all the derived quantities. Table [3] gathers the median error 
for all the sources. It is seen that the typical errors are very 
small. This is caused by the fact that the Kron apertures de- 
ployed by SEXTRACTOR include of the order of lO^pixels, 
and hardly ever less than 10^ pixels. The expected errors in 
the structural parameters are therefore very small. The set 
of table entries above the horizontal line in Table [3] are the 
statistical uncertainties. Table entries below the horizontal 
line gather the systematic uncertainties for the asymmetry 
indices, derived by calculating the asymmetry indices for the 
GALFIT models. 



4 MERGER IDENTIFICATION USING 
STRUCTURAL PARAMETERS OF THE 
RESIDUALS 

As it was mentioned in JT] the merger fraction is defined 
as the fraction of galaxies with an ongoing merger episode 
that is found in any given sample, which is often selected as 
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a mass- limited one. This is a fundamental issue to galaxy 
evolution studies and therefore the automated identification 
of mergers is of great importance. This section presents the 
main contribution from this paper. It shows that the use of 
the structural parameters of the residual images allows to 
identify mergers like the use of the structural parameters of 
the galaxies themselves. It also shows that the merger sam- 
ples obtained using the properties of the residual images are 
of better or comparable statistical quality than the samples 
that are culled using the morphological information of the 
original images. The precise definition of the statistical qual- 
ity of a sample, which is given in ^4.11 is used together with 
the set of galaxies defined in ^2.11 that were morphologically 
classified by the STAGES team. This key ingredient allows 
us to obtain a sample of visually detected mergers, which 
will be used as a training set for the method presented in 
this section. Many studies have used the morphological or 
structural properties of galaxies to estimate the merger frac- 
tion. Usually, these works make use of the CAS system or of 
the G-M20 system. Eac h of these systems has its own pros 
and cons. For instance, IConsel ice. Raigor. Mvers' (2008|) 
concludes that the G-M20 system discovers more mergers 
than the CAS methodology, although it also picks up more 
interlopers. It is however the case that all the structural ap- 
proaches that have been presented so far have only taken 
advantage of the structural parameters of the real, direct 
images of galaxies. In addition, CAS t ypically re covers a 
fairly hi^h fract ion (50% to 70% in iJogee et al.1 (|2009l ): 
iHeiderman et al.l ([2009, \) of visually classified mergers, but 
it is also significantly contaminated by dusty, highly inclined 
non-interacting galaxies. These latter galaxies have low level 
as ymmetries caused b y star formation episodes, as was noted 



as ymmetries causea p y 
in bogee et all (|2009l ). 



4.1 Statistical Quality of Samples 

In science, one often confronts the problem of finding an al- 
gorithm or method to select a sample of items from a larger 
parent population with the condition that the selected items 
have to satisfy some requirements of scientific interest. How- 
ever, one rarely has a mechanism to retrieve all the items in 
the parent population that satisfy the needed requirements, 
and it is also very unlikely that the method is able to retrieve 
the required items only. One is then forced to speak about 
the sensitivity and specificity of the selection process. 
The sensitivity, also known as the recall ratio, is defined 

as: 



#True Positives 



#True Positives + # False Negatives 



(11) 



This is more commonly known as the completeness in astro- 
nomical literature. 

The specificity is defined as: 



P ■ 



#True Negatives 



#True Negatives + # False Positives 



(12) 



In the above definitions, a "True Positive" is a recovered 
item that did indeed present the required properties. A 
"False Negative" is an item that was not retrieved by the 
culling algorithm but did present the needed properties. 
These latter errors usually refiect an excessive skepticism. 
A "True Negative" is an item that was rightfully rejected 



by the selecting process since it did not have the required 
properties. A "False Positive" is an item that was incorrectly 
picked up by the sampling algorithm, but that does not have 
the properties of interest. 

The sensitivity and the specificity can b e combined into 
a single number, known as the F — score, F/s (|van RiisbergenI 
1979). This is a measure of sample purity, and it is just a 
weighted harmonic average of r and p. 



(1 + X p X r 
X p + r) 



(13) 



where /3 is a control parameter that regulates the relative 
importance of r with respect to p. This is a user-supplied 
value that depends on the particular goals of the test|j. In 
this work, a value of /3 = 1.25 is used, which can be thought 
of as weighing completeness more than the lack of contam- 
ination. The use of this value will be justified in ^4.41 This 
choice leads to a galaxy sample that contains most mergers 
from that training set, although the corresponding contam- 
ination is rather high. This will be further discussed in ^5.11 

The F — score is used in the current study in order 
to grade the performance of a number merger diagnostics 
at separating a merger sample from its parent population. 
Galaxies undergoing a merger episode play the role of the 
"items presenting the required properties" discussed above, 
and the parent population used here is of course the galaxy 
sample defined in ^2.11 

In this context, a "merger diagnostic" is defined as a 
two-dimensional diagram in which the parent population of 
galaxies described in ^2.11 is presented. In these plots, the 
horizontal axis is one structural parameter and the vertical 
axis is another morphological parameter, both selected from 
the set of indices described in ^3.31 In these diagrams, merger 
galaxies should preferentially occupy specific regions. For in- 
stance, mergers should have large asymmetries and higher 
than average values of G. This is exploited by searching for 
the best border that separates mergers from other galaxies 
in each of these diagrams. The "border" of a diagnostic is 
defined as a second order polynomial in the horizontal co- 
ordinate that maps to the vertical coordinat^fl- Although 
it would be possible to use this same technique with more 
than two parameters, this would make it more difficult to 
interpret the resulting 3-D parameter spaces. For this rea- 
son, the merger diagnostics considered in the current work 
are simply 2-D. 

Galaxies are then classified into four different types, de- 
pending on the side of the border in which they fall and on 
whether they are involved in a merger or not. 

• Mergers that fall in the merger side of the border are 
the "True Positives". 

• Mergers that do not fall in the merger side of the border 
are called "False Negatives". 



^ The reader can probably work out what is the value of ^ used 
by the managers of airport security screenings, where false pos- 
itives represent a small additional test but false negatives have 
disastrous consequences. 

^ The Gini-M2o method, in which galaxies with too high values of 
G for their M20 value are classified as mergers, is a basic example 
of this approach. 
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Parameter 1 . 



Figure 1. Schematic graphical explanation of the F — score 
method used. The horizontal axis is one of the structural param- 
eters calculated either for the original images or for the residual 
images. The vertical axis is another morphological parameter, re- 
gardless of its nature. The larger dots represent the galaxies that 
were marked as mergers in the STAGES morphological catalog, 
and the smaller dots are galaxies that were not classified as merg- 
ers. The thick line is then the border, which is fine tuned by the 
Amoeba algorithm so that most of the large dots lie above it, 
while at the same time most of the small points are located un- 
derneath it. In this representation, the merger side of the border 
is clearly the zone of the diagram above the border. 



• Non mergers that however fall in the merger side of the 
border are regarded as "False Positives" . 

• Non mergers that do not fall in the merger side of the 
border are of course "True Negatives". 

In all the above definitions, the merger side of the border is 
to be understood as the zone in the diagram in which the 
majority of mergers exist. 

The best border is then the border that maximizes the 
F — score. In this step, mergers serve as buoys or perhaps 
better as a training set in order to find the best border. This 
maximizatio n is done by mean s of the Amoeba algorithm 
explained in IPress et al.l (| 19881 ) , using the polynomial coef- 
ficients as the problem parameters. 

The method is schematically presented in Figure [T] 
Both axes represent a dummy structural parameter, which 
could be any one of the structural parameters defined in 
^3.31 The large dots represent merging systems and the small 
dots are non mergers. The thick line is a best border found 
by the F — score maximization algorithm. This is a second 
order polynomial in the Parameter #1 into the Parameter 
#2 dimension. 

^4.31 presents this optimization for a number of merger 
diagnostics. Some of the merger diagnostics presented are 
taken from the literature, while others are introduced there 
for the first time. The optimization process used tries to 
maximize the F — score number, which is derived from the 
r and p statistics. A meaningful assessment of the resulting 
contamination requires an accurate estimate of the merger 
fraction in the galaxy sample used here. The contamination 
ratios for the best performing diagnostics is presented in 



4.2 The Training Set for the F — score Technique 

The next question is then what objects from the parent pop- 
ulation of galaxies are to be taken as true mergers and there- 
fore used as the training set for the F — score maximization 
technique, presented in ^4.31 To this end, the morpholog- 
ical information presented in the STAGES morphological 
catalogue is used. 

For the purposes of this study, galaxies classified as 
mergers by at least two of the visual observers of the 
STAGES team are defined to be as mergers. This subset 
includes 39 objects and constitutes the training set. On the 
other hand, sources which were classified as mergers by one 
or none of the STAGES observers were then considered as 
non-mergers. This subset is made of 1498 sources. This sub- 
set includes 83 sources that were regarded as mergers by 
only one of the STAGES visual classifiers. These objects 
were not included in the training set in order to obtain a 
more robust merger training sample. 

Figure [2] is an image atlas showing images of all the 
galaxies from this training set. This figure contains both 
the angular and physical scale of each inset, the STAGES 
ID of each source, the environment (Cluster "C" , and Field 
"F"), and the morphological type (Elliptical "E", Lenticular 
"SO", Spiral "Sp", and "0th"). Each inset also includes the 
number of observers that agreed on whether that particular 
galaxy is a merger or not. 

Figure [2] shows that most of the 39 galaxies used as 
the training set are indeed bona- fide mergers. It is also true 
that a small fraction might be more marginal. These latter 
galaxies need to be included in the merger set because one 
of the main goals of the current work is to explore whether 
the analysis of the structural parameters of the residuals can 
be used to detect minor mergers reliably. 

4.3 Comparison Between Merger Diagnostics. 

Given the wealth of structural parameters that have been 
calculated for the galaxies in the parent population, it is 
better to start presenting how the F — score method works 
for the G'(Obj) — M2o(0bj) plane. This is a very well known 
merger diagnostic that has already been used in the litera- 
ture. In this subsection, the image for which the structural 
parameters are calculated is denoted inside parentheses af- 
ter the name of the structural parameter itself. "Obj" means 
that the morphological parameter was obtained in the orig- 
inal image, "MDL" refers to the Sersic model, and "Res" 
implies that the parameter was obtained for the residual 
image. 

Figure [3] shows the G'(Obj) — M2o(0bj) plane used as a 
merger diagnostic. Large symbols are the galaxies that were 
marked as mergers by the STAGES observer^, and the 
smaller symbols show the galaxies that were not regarded 
as mergers. Cyan filled circles represent irregular objects, 
blue squares give the location of spiral galaxies, black trian- 
gles represent lenticular systems and red diamonds denote 
elliptical galaxies. This panel shows the "best" border as 
a thick, green solid line. The dashed black line is just the 



^ This is thus the training set the Amoeba algorithm will use to 
find the best border to separate mergers from non-mergers. 
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Figure 2. Image atlas of the merger training set. The COMBO-17 ID is included in each panel. This figure shows the objects which 
were classified as mergers by at least two of the STAGES team visual observers. Each panel shows three different insets. The first image 
is the direct image, the second image is the model created by GALFIT, and the third image is the residual image. This latter panel is 
shown with an inverted look-up-table whose dynamic range is 15% that of the other two images. This is don in order 
to enhance the visibility of the fainter features in the residual images. 
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Figure 2 continued Atlas of the training set galaxies. 
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Figure 2 - continued Atlas of the training set galaxies. 





initial border that the Amoeba algorithm is given to start 
its iterations. This is just a rough guess given the location 
of the larger symbols in the diagram. The final best bor- 
der the Amoeba algorithm obtains does not depend on the 
initial guess as long as this initial guess is reasonable. The 
polynomial that defines the best border is also given within 
the figure, together with the resulting completeness, con- 
tamination, and F — score values. Here, the completeness 
is defined as the number of clear visual mergers above the 
border, divided by the total number of visual mergers (i.e., 
the sensitivity). 

It is seen that the F — score maximization algorithm has 
been able to find most of the merging systems that consti- 
tute the training set. The statistical quality of this sample 
is F = 0.77. The completeness is r = 0.79, and the contam- 
ination by objects not classified as mergers is 1 — p = 0.28. 
It is interesting to note that this method rejects most of the 
lenticular objects from the training sample. 

The next merger diagnostic presented is the A(Obj) — 
RFF plane. This merger indicator is motivated in the more 
common A(Obj) — 5'(0bj), using the RFF instead of the 
dumpiness since these two quant ites are very similar. Fig- 
ure 2] presents this test. This is the first test that makes use 
of the structural parameters of the residual images. In this 
test, the initial border is A(Obj) = 0.30, simply because the 
asymmetry is expected to bear the highest predictive power 
in this test. 

Figure m shows that the statistical quality of the sample 
obtained using this criterion has improved significantly with 
respect to the results achieved by the G'(Obj) — M2o(0bj) 
diagnostic. The sample purity is F = 0.85, the sensitivity 
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Figure 3. The well-known G(Obj) — M2o(Obj) merger criteria 
applied to the STAGES sample. Large symbols are sources that 
were marked as mergers by the STAGES team observers. Blue 
squares show spiral galaxies, black triangles denote lenticular sys- 
tems and red diamonds present elliptical galaxies. Irregular and 
disturbed systems are represented as beige filled circles. The green 
line is the best border found by the Amoeba algorithm, and the 
black, dashed line is the initial guess this algorithm is given. The 
optimal value of F — score, together with the sample completeness 
and specificity are given within the figure. 

is r = 0.92, and the specificity is p = 0.76. Also, Figured 
shows a clear correlation between A(Obj) and the RFF. It 
might be obj ected that the u sual limit in the asymmetry 
introduced in IConselicl (|2003l ) is A(Obj) > 0.35, while the 
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Figure 4. The A(Obj) — RFF merger diagnostic. Symbols and 
information as in figure [3] 
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Figure 5. This figure presents the G(Obj) — A(Obj) plane used 
as a merger diagnostic. Symbols and information as in figure [3] 



limit suggested by this test is A(Obj) > 0.20. This is not 
only caused by the training sample used, but by the current 
choice of the /3 parameter, which is designed to weight com- 
pleteness more than specificity. Had these condition been dif- 
ferent, the resulting best borders would have changed. This 
method therefore imperatively requires an objective calibra- 
tion in order to produce meaningful, physically motivated 
borders and hence reliable merger fractions. This objective 
calibration will be produced in a forthcoming paper, using 
full fledged N-body simulations of galaxy mergers. Also, it 
has to be borne in mind that the A(Obj) > 0.35 criterion 
is tuned to detect major mergers, and one of the aims of 
this work is to improve the morphological detection of mi- 
nor merger episodes. Figure 2] also shows that the F — score 
maximization algorithm has indeed found that the "best" 
border is very different from the initial, fiat guess. In par- 
ticular, the Amoeba algorithm has discovered the correla- 
tion between A(Obj) and RFF^ and takes advantage of it 
by converging towards a line that cuts the correlation in a 
perpendicular way. It is also tempting to think that the use 
of B+D decompositions could produce better overall fits for 
the SO and Sp galaxies, leading to lower RFF numbers. This 
in turn could reduce the number of false positives for these 
sources, increasing the potential of this structural parameter 
as a merger diagnostic. This will be explored in the future. 

The results shown in Figures [3] and S] make it possible 
to think that a good merger diagnostic could be put to the 
test by combining G'(Obj) with A(Obj). Figure [5] presents 
this investigation, which confirms the previous ideas with a 
sample purity oi F — 0.86, a recall r = 0.90 and a fairly high 
specificity p — 0.82. The algorithm has correctly identified 
that the best border in this plane is a diagonal, as expected. 
The final border is found to be very close to initial border 
tried by the minimization algorithm. 

The following merger criterion considered continues 
with the exploration of the morphological parameters of 
the residual images. Figure [6] shows the G'(Res) — M2o(Res) 
merger test plane. This diagnostic was motivated by two 
main ideas. 

• If the galaxy that was fitted and removed by the smooth 
Sersic model was indeed involved in a merger episode in a 
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Figure 6. The G(Res) — M2o(Res) plane as a merger diagnostic. 
Symbols as in Fig. [3] 

very late stage, the residual image should expose the effect 
of the fainter component as a bump. This bump will then 
be easily detected in G'(Res) because it will form a reduced 
number of high- intensity pixels, surrounded by a large num- 
ber of pixels with very low intensity values with an aver- 
age value of 0.0. This will boost the value of G'(Res). If, on 
the other hand, the galaxy is well described by the smooth 
model, the pixel intensity of the residuals will all cluster 
around 0.0, and G'(Res) will be very close to 0.5. 

• In the case of the horizontal axis, which is M20 (Res) , the 
situation is very similar. If there is an off-centre bump in the 
residual image, it will probably be among the 20% brightest 
pixels of the residual image. This will enhance the value of 
M20 with respect to the situation in which no substructure 
is seen in the residuals. In this latter situation, the brightest 
20% pixels of the residuals will be randomly distributed in 
the image, albeit with a preference for the central values. 

Figure [6] confirms all these expectations. It is seen that 
the sample purity is F = 0.79, with r = 0.77. The specificity 
is p = 0.82. It is also clearly seen that the predicting power of 
this diagnostic is mostly associated with G'(Res). The results 
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Figure 7. The G(Res) — A(Obj) plane as a merger diagnostic. 
Symbols as in Fig. [S] 



Figure 8. The RFF — G(Res) plane, used as a merger criterion. 
Symbols as in Fig. [3] 



of this test suggest the use of the G'(Res) — A(Obj) plane as 
a merger diagnostic. This idea is presented in Figure [71 The 
correlation between the two quantities presented is evident 
in this figure. This alone indicates that G'(Res) is bound to 
be a good merger tracer, even when used by itself. For this 
test, the sample purity is F = 0.85, with a good complete- 
ness r = 0.90 and a better specificity p = 0.78. This merger 
indicator is then comparable to the G'(Obj) with A(Obj) 
indicator. 

The next natural merger criterion that the conclu- 
sions drawn from Figures [6] and U lead to explore is the 
G'(Res) — RFF plane, shown in Figure (8] These two param- 
eters have been selected because the RFF correlates with 
A(Obj). It is therefore expected that this indicator will also 
be able to separate mergers from other galaxies. The G(Res) 
was used because Figure [6] shows it could single out merg- 
ers almost as a standalone indicator. Figure [8] shows that 
the RFF is indeed correlated with G'(Res), and that this 
correlation could be exploited to identify mergers. The sam- 
ple quality is F = 0.84, with a specificity p = 0.80. The 
completeness is r = 0.87. It is however the case that the 
best border in Figure [8] is approximately horizontal, indicat- 
ing that the Amoeba algorithm does not take advantage of 
the correlation between the RFF and the Gini index of the 
residuals. 

This work would not be complete without studying 
the CAS parameter space. Figure [9] explores the A(Obj) — 
5'(MDL), A(Obj) - C(MDL), and the A(Obj) - C(Obj) 
planes. In all the panels in Figure (9] the vertical axis is 
A(Obj). However, the horizontal axis of the upper panel 
is the dumpiness 5'(0bj), the horizontal axis of the mid- 
dle panel is C(MDL) and the horizontal axis of the lower 
panel is C(Obj). In all these cases, the initial guess line 
is obviously motivated by the classical A(Obj) > 0.35, 
A(Obj) > 5'(0bj). It is seen that the classical CAS crite- 
rion works very well, achieving a sample purity F = 0.86, 
with a good recall r = 0.95 and a fairly high specificity 
p — 0.76. It is also seen that the Amoeba algorithm has 
found that the clumpier mergers need also be more asym- 
metric in order to be classified as such. Comparison between 
the A(Obj) -C(MDL) and A(Obj) - C(Obj) highlights that 



the use of a smooth model to calculate C(MDL) is slightly 
beneficial for merger detection. This comes from the higher 
sample purity obtained with the A(Obj) — C(MDL) plane. 
For this latter test, F = 0.82, with a specificity p = 0.68. 
The completeness is r = 0.95. This is clearly a success of 
the A(Obj) parameter. The contamination rate obtained by 
the use of the A(Obj) — 5'(0bj), as presented here, will be 
shown in ^5.11 

It is interesting to complete the analysis of the CAS pa- 
rameter space by deriving the values of the r, p, and F, pa- 
rameters for the A(Obj) ^ 0.35, A(Obj) > 5'(0bj) criterion, 
which is a frequently used criterion for CAS based major 
merger studies. Using these limits, r = 0.21, p = 0.98, and 
= 0.30. It is clear that with this criterion the merger sam- 
ple obtained discards the majority of the non-mergers, but 
is very incomplete. The contamination of this latter sample 
will be presented in ^5.1[ but it has to be borne in mind that 
the optimization used in this work aims towards complete- 
ness, while the goal of the A(Obj) ^ 0.35, A(Obj) > 5'(0bj) 
limits is a clean sample of major mergers. The two diagnos- 
tics cannot be compared directly, then. 

The results presented so far suggest that the asymme- 
try is a very good merger indicator. It is also the case that 
the G'(Res) and the RFF can also be used as merger indica- 
tors by themselves. In the spirit of this work, the following 
planes studied will include the asymmetry of the residuals 
A (Res). This morphological parameter is explored because 
the asymmetry of a system is expected to be boosted after 
the subtraction of an intrinsically symmetric profile such as 
the Sersic profile used here. Figure [TOl shows the success of 
this approach, in particular on its lower panel. 

Both diagnostics shown in Figure [TOl are found to work 
very well. In particular, the RFF — A{Res) produces the best 
sample purity of all the diagnostics tested here that exploit 
the structural parameters of the residuals, with a very high 
completeness. As expected, the asymmetry of the residuals 
is boosted with respect to the asymmetry of the original im- 
ages. The best border for the RFF — A(Res) diagnostic is no 
longer horizontal, implying that A (Res) has enough predict- 
ing power so as to curve the best line. The best border line 
of the RFF — A (Res) test appears to be more curved than 
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Figure 9. Upper panel: The A(Obj) — 5'(Obj) plane as a merger 
diagnostic. Middle panel: The A{Ohi) — C(Mdl) indicator. Lower 
panel: The A(Obj) — C(Obj) merger test. Symbols as in Fig. |3] 

that of the G(Res) — ^(Res) merger test. This is consistent 
with the previous conclusion that G(Res) could also be used 
as a standalone merger diagnostic. The purity of the sample 
obtained using the G(Res) — ^(Res) indicator has F = 0.84, 
with a recall r = 0.90 and a specificity p = 0.78. The purity 
of the sample obtained using the RFF — j4(Res) indicator 
has f = 0.86, with a r = 0.92 and p = 0.78. This is then 
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Figure 10. Upper panel: The A(Res) — G(Res) merger test. Lower 
panel: The A(Res) — RFF plane as a merger diagnostic. Symbols 
as in Fig. [3] 



the best test, and this confirms that the structural param- 
eters of the residuals can indeed be competitive if used as 
merger diagnostics. The corresponding contamination stem- 
ming from the RFF — A (Res) diagnostic will be shown in 
^5.1[ where a sample of visually classified mergers will be 
used to establish the merger prevalence. Whether the use of 
these parameters can be used to probe deeper in the lumi- 
nosity function in order to detect mergers with mass ratios 
larger than 10:1 or if these parameters allow to trace the 
merger event up to later stages in which the less luminous 
galaxy is almost engulfed by the host will be presented in a 
forthcoming paper. 

The quality of the merger samples obtained using the 
RFF — A (Res) test is thus the best one of all the residual- 
based merger diagnostics that have been explored in this 
work. The discussion presented in JS] will be based on the 
merger sample obtained using this merger test. Also, ^5.21 
shows some objects that are thought to be involved in minor 
merger events that have been detected in the RFF — A(Res) 
plane as mergers that would have been missed by the usual 
CAS criterion A(Obj) ^ 0.35, A(Obj) > 5'(0bj). 

Table |4] gathers a summary of the different merger di- 
agnostics used and their respective statistical performances, 
for convenience. 
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Figure. 


F — score 
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G(Obj)-M2o(Obj) 


El 


0.77 


0.79 


0.82 


A(Obj) - RFF 
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0.85 


0.92 


0.76 


G(Obj) - A(Obj) 
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0.86 


0.90 


0.82 


G(Res) - M2o(Res) 
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0.79 


0.77 


0.82 


G(Res) - A(Obj) 


m 


0.85 


0.90 


0.78 


RFF - G(Res) 


m 


0.84 


0.87 


0.80 


A(Obj) - ^(Obj) 


m 


0.86 


0.95 


0.76 


A(Obj) - C(Mdl) 


m 


0.82 


0.95 


0.68 


A(Obj) - C(Obj) 


m 


0.77 


0.72 


0.86 


A(Res) - G(Res) 


Co] 


0.85 


0.90 


0.78 


A(Res) - RFF 


Co] 


0.86 


0.92 


0.78 



Table 4. Summary of the different merger diagnostics tried and 
their statistical performance. 



4.4 The Effect of the /3 Parameter. 

The F — score maximization technique used in the current 
work is simply a way to select a number of objects from a 
parent population whose structural properties are similar to 
those of the training set. The method is thus only as good 
as the training set. The /S parameter in the definition of F^ 
is the ingredient used by this statistic to decide whether or 
not the structural properties of an object are close enough 
to the structural properties of the objects in the training 
set to be considered as a merger candidate. The /3 param- 
eter is therefore determined by the scientific needs of the 
sampling process. Figure CD presents a recalculation of the 
"best" border of the ^(Res) — RFF diagnostic using two 
different values of /3, /3 = 2.0 (upper panel) and /3 = 0.5 
(lower panel). 

Comparison of Figure CD with Figure CHI clearly indi- 
cates that the /3 parameter has a decisive impact on the sam- 
ple of potential mergers obtained by the Amoeba algorithm. 
If completeness is considered to be much more important 
than specificity, as is done in the upper panel of Figure CD 
the Amoeba algorithm increases the recall ratio by 3% with 
respect to the value achieved in Figure CHI This amounts to 
the inclusion of just an additional object from the training 
sample, at the price of littering the sample with 150 non- 
mergers. This decreases the precision by 10%. On the other 
hand, the lower panel of Figure CD shows what happens if 
the specificity is weighted more than the completeness. It is 
first seen that the new merger diagnostic excludes a higher 
number of non-mergers, raising the precision by 12%. This 
comes at the expense of missing 11 objects from the training 
set. The choice of (3 = 1.25, adopted in ^4.11 is thus a good 
compromise between these two options. 

It is highlighted here that no single value of /3 can be 
considered to be "correct" . The value of this parameter is set 
by the goals of the test. If, for instance, the objects selected 
by the method are to be the targets of a spectroscopic follow- 
up programme where targeting only true mergers is deemed 
essential, it would be advisable to give a higher weight to 
the specificity. 
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Figure 11. Recalculation of the A(Res) — RFF merger diagnostic 
using /3 = 2.0 (upper panel) and f3 = 0.5 (lower panel). Symbols 
as in Fig. [3] 



5 VISUAL ASSESSMENT OF THE MERGER 
SAMPLES OBTAINED. 

As it is clear from the previous discussion, the method pre- 
sented in the current work is optimized to be very complete, 
hopefully detecting minor mergers thanks to the use of the 
structural parameters of the residual images. It is therefore 
needed to establish both the contamination by non-mergers 
and the success in recovering minor mergers. The study of 
the contamination by non-mergers is gathered in ^5.11 On 
the other hand, ^5.2| presents a number of examples of minor 
merger candidates that have been detected through the use 
of the structural parameters of the residuals that would oth- 
erwise have been missed by the traditional CAS diagnostics. 



5.1 Contamination by Non-Mergers. 

The different galaxy samples obtained by the sheer use of 
the various indicators that have been put to the test in ^4.31 
need to be further evaluated. It is clear that not all the ob- 
jects that test positive (i.e., sources that fell in the merger 
side of the merger tests used) to these criteria can be merg- 
ers. In particular, the contamination is the key statistic that 
needs to be included in this appraisal. This section presents 
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the non-merger contamination of the galaxy sets that are ob- 
tained from the use of the merger diagnostics that attained 
the highest F — score values in ^4.31 

The contamination is here defined as: 

^ _ #Non — mergers that however test positive. 
~ #A11 positives. 

where the objects that test positive are those in the merger 
side of the "best" borders calculated above. The denomina- 
tor of this fraction includes both mergers and non-mergers, 
which implies that the contamination depends on the merger 
fraction in a non-linear way. 

The next step is thus oriented towards obtaining a com- 
plete and accurate estimate of the fraction of mergers found 
in the parent population studied. The total number of good 
mergers is here calculated by adding together the objects 
included in the training set of galaxies used above, and a 
number of merger galaxies that were recovered during an 
additional visual classification which will be described be- 
low. This further observational classification is justified be- 
cause the STAGES morphological catalogue used to define 
the training set of objects for the F — score maximization 
technique is a general morphological catalogue that is not 
designed to split the parent population into mergers and 
non-mergers. In particular, the STAGES observers were not 
specifically looking for the minor mergers whose detection 
is the goal of this study. This additional study thus serves to 
check whether minor mergers do share the structural proper- 
ties of major mergers, which is the main assumption behind 
the use of the F — score number as a diagnostic discrimina- 
tor. 

The new visual assessment examines the set of galaxies 
obtained by the blind application of the A (Res) — RFF di- 
agnostic as shown in Figure [101 This merger test was shown 
in ^4.3l to yield the highest F — score and specificity num- 
bers of all the diagnostics that make use of the structural pa- 
rameters of the residual images and is therefore more likely 
to produce a clean list of mergers. This set of galaxies is 
made of 36 = 0.92 x 39 objects from the original training 
set and 332 = (1 — 0.78) x 1498 objects that were not in- 
cluded in the training set. These latter systems are the False 
Positives involved in the calculation of the p statistic. The 
majority of these objects (282) have zero merger marks in 
the STAGES morphological catalogue and 50 of them only 
have one merger mark. These simple statistics motivated us 
to examine more closely the 332 objects not included in the 
original training set that however fell in the merger side of 
the "best" border in the A (Res) - RFF plane. The pur- 
pose of this further investigation is to establish whether or 
not those objects could be mergers that escaped the original 
assessment of the STAGES team observers. This will allow 
an accurate and non-parametric determination of the merger 
prevalence in the parent population of galaxies studied. To 
this end, four of the authors of this paper (A AS, CH, EFB, 
MEG) reinspected the 332 False Positives together with 332 
randomly selected True Negatives as a control sample. These 
four independent assessments were then combined into a sin- 
gle trinary outcome, splitting these 664 sources into three 
different sets, (i) clear mergers, (ii) clear non-mergers, and 
(iii) the dubious cases. These 664 objects were then placed 
in the A(Res) — RFF plane. This is presented in Figure [T2l 

Figure [12] immediately shows a correlation between the 




Figure 12. Visual assessment of the galaxy sample obtained us- 
ing the RFF — A(Res) merger test. Red, solid dots represent the 
location of the new best merger cases (55 points), blue diamonds 
give the location of the dubious cases (41 points) and the open, 
black points show the location of the non mergers (568 points). 
The black line shown is the same line shown in Figure [To] 

location of the mergers that have been recovered by the new 
visual merger re- assessment and the RFF — A (Res) "best" 
border. It is seen that the majority of the recovered mergers 
and dubious mergers lies above the diagnostic line. It is also 
highlighted that none of the objects shown in Figure [T2l 
belongs to the original training set. 

Figure [13] presents eight of the galaxies that were 
reinspected during this further classification. This figure 
presents four rows of two objects each, marked with their 
respective COMBO- 17 IDs. The upper row shows the ob- 
jects that tested as positives to the RFF — A{Res) diagnostic 
and were subsequently classified as mergers in the new vi- 
sual assessment. The second row presents mergers that were 
regarded as such that however test negative to the diagnos- 
tic. These latter objects are very rare. The third row is made 
of non-mergers that nevertheless tested positive in the au- 
tomated test, and the bottom row presents the non mergers 
that fell below the "best" border in Figure [TOl The percent- 
ages given in each row are the frequencies with which each 
of the different possibilities appears. These percentages do 
not add up to 100% because the galaxies that were classified 
as dubious mergers are not shown. 

The new visual inspection has recovered a total of 55 
clear merger systems together with 41 dubious cases. The 
remaining 568 objects were classified as non mergers. It is 
interesting to note that 36% of the 55 clear merger cases 
had just one merger mark in the original STAGES mor- 
phological catalogue, 18% of the new possible mergers were 
considered as mergers by only one of the STAGES classi- 
fiers, and only 6% of the non-mergers had received a single 
merger mark. This adds up to a total of 60 objects with 
only one merger mark in the whole pool of the 664 objects 
inspected. A total of 50 sources out of these 60 objects are 
found above the "best" border defined in Figure [TO] for the 
RFF - A{Res) merger test. In addition, the RFF - A{Res) 
merger test as presented in Figure [10] detects 93% of the re- 
covered clear merger cases, 88% of the new dubious mergers 
and 43% of the non-mergers. 
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Figure 13. Examples of galaxies examined in the new visual 
assessment. COMBO-17 IDs are shown for each galaxy. The per- 
centages given reflect the frequencies of the four different possi- 
bilities shown. 



Thus, the total number of objects selected by the 
RFF — A (Res) test can be broken down in the following 
way: 

(i) The 36 = 0.92 x 39 objects included from the training 
set. 

(ii) The 51 = 0.93 x 55 new mergers recovered by the new 
visual inspection. 

(iii) An unknown number of merger systems between 
and 36 = 0.88 x 41 that were classified as dubious mergers 
by the visual re- assessment. Under the assumption of a flat 



probability distribution this number could be represented as 
18 ±10. 

(iv) A total of approximately 260 = 332 — 57 — 19 non 
mergers. 

Therefore, the total number of bona- fide mergers above 
the RFF — A (Res) "best" border line is 105, and the total 
number of objects of all classes above this line is 368 = 332 + 
36. The final contamination ratio is 71% = (368— 105)/368, 
which is fairly high since it means that 70% of all the objects 
set aside by the blind use of the diagnostic are non-merger 
contaminants. This number is close, but conceptually differ- 
ent from the fraction of non-mergers that are found to test 
positive in during the additional visual classification, which 
is 74%. 

In the same way, the non-merger contamination associ- 
ated with the A(Obj) — 5'(0bj) test as shown in Figure [9] is 
71±3% = (397-(37+(360/286)x(49+(0.5±0.3)x28)))/397. 
This calculation takes into account that there are 397 = 
39 X 0.95 + 1498 x (1 - 0.76) objects above the "best" border 
defined in Figure [9] and that not all the dubious mergers will 
indeed be mergers. This latter consideration is again made 
assuming a flat probability distribution. This contamination 
that affects the sample selected by this method is seen to be 
fully compatible with the contamination that is calculated 
for the RFF — A (Res) test. If a similar analysis is carried 
out for the A(Obj) — 5'(0bj) criterion using the traditional 
CAS limits A(Obj) ^ 0.35, A(Obj) ^ 5'(0bj), the contami- 
nation is 50% = 19/(19 + 17). Thus, the contamination by 
non-mergers found in this set of galaxies is lower than the 
one found for the other two merger diagnostics explored, but 
its completeness is obviously much lower. 

The next step in this analysis is then the study of the 
negative detections, which focuses on the 1169 sources be- 
low the "best" border line presented in Figure [10] These 
1169 galaxies include 3 galaxies from the original training 
set (these are the "False Negatives" in the F — score analy- 
sis). These 1169 galaxies will also include an indeterminate 
number of mergers that can be estimated from the visual 
re-assessment results by multiplying the fraction of merg- 
ers found in the 332 objects from the control sample by 
the total number of objects below the "best" line (1169). 
The new visual classification discovered 4 clear merger and 
5 dubious cases in the 332 objects in the control set. If 
these latter sources are given a weight of 0.5, the fraction of 
mergers below the "best" line of the RFF — A (Res) test is 
2.2 ± 0.4% = (3 + (1169/332) x (4 + (0.5 ± 0.3) x 5))/1169, 
where the error interval again assumes a flat probability 
distribution for the dubious mergers. If the galaxies from 
the original training set and the newly identified dubious 
mergers are not included in this calculation, the fraction 
drops to 1% = ((1169/332) x 4)/1169, which is the num- 
ber given in Figure [13] for the negative contamination ratio 
of the visual re- assessment. The negative contamination ra- 
tio is then very low, indicating that this technique is very 
powerful as a negative merger test. Furthermore, if this neg- 
ative contamination ratio is derived for the A(Obj) — 5'(0bj) 
test as presented in Figure [9] the percentage is 4.0 =b 1% = 
(2+(1140/378) X (6+(0.5=b0.3) x 13)) /1140, which is compat- 
ible but slightly worse than the result for the RFF — A{Res) 
test. This latter calculation uses that there are 378 objects 
in the total pool of objects inspected during the second vi- 
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sual classification that fell below the "best" border shown 
in the upper panel of Figure [9l Finally, the corresponding 
negative contamination for the traditional CAS criterion is 
5% = 77/(1424 + 77). This is compatible with the negative 
contamination for this diagnostic as presented in Figure [9l 
In summary, despite the fact that the contamination ra- 
tio is fairly high for the positive detections, it is very low for 
the negative detections. The above considerations therefore 
lead us to conclude that the RFF — A (Res) minor merger 
diagnostic presented in Figure [10] works best as a negative 
test. In particular, it could be possible to use this auto- 
mated technique wit h large area surveys such as th e APM 
Maddox et atlll QQOV 2dFGRS (iFolkes et al.lll999l l. SDSS 



Abazaii an et al. 2009). UKIDSS (IPve et al.ll2Q06l), KID^B, 



VIKI NCQ, GAM A (IPriver et al.ll2009l), GEMS (|Rix et al.l 
l2004l ). CO SMOS (Ik oekemoer k Scovilld [iooHl l.DE^. or 
the LSST (llvezic et a l. 2008) by calibrating the "best" bor- 
ders to the new observational conditions and wavebands us- 
ing a reduced and manageable number of objects within a 
predefined redshift window up to a certain magnitude. These 
new "best" borders should then be used to split the tar- 
get galaxy population. The technique illustrated here would 
then produce two different sets of galaxies. One set would be 
almost completely free of mergers and the other set would 
include the overwhelming majority of mergers. This latter 
sample would need further purification in order to produce 
a clean sample of mergers. This recalibration step is particu- 
larly needed in the case of surveys including U band observa- 
tions. In this case, even very minor wet mergers will leave a 
larger impact on the general appearance of their host galax- 
ies because the youngest stars will be much clearly seen. 
This step could also help understand the effect of the pho- 
tometric band chosen on the optimization process used here. 
It is clear that the additional pruning can not be done us- 
ing the structural properties of the galaxies only, since these 
have been fully exploited here. In addition, one clear conclu- 
sion from this visual assessment which will be strengthened 
in ^5.21 is that even a visual inspection cannot unambigu- 
ously tell whether a particular object is involved in a merger 
episode or is just the product of a by-chance alignment or an 
HII region. On the contrary, additional information such as 
colours or kinematical information obtained using Integral 
Field Units spectrographs needs to be included. Note that 
this problem affects all automated methods based on struc- 
tural parameters that we have explored in this paper and 
are commonly used in the literature. This issue is merely a 
consequence of the empirical existence of visually classified 
mergers with undisturbed r norph ologies, which was also ob- 
served in Heiderman et al.l (|200 9). The above considerations 
make it advisable to seek an inclusive criterion, with a high 
recall ratio, rather than a specific one, which is refiected in 
the choice of the /3 parameter. 



5.2 The Search for Minor Mergers. 

As it has been mentioned, the main driver behind the very 
high completeness that the method presented here is tuned 
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to achieve is the detection of minor mergers. It is therefore 
interesting to test the ability of the methodology introduced 
here to detect minor mergers. This is done by selecting a 
set of visually identified mergers that would not have been 
detected by the usual cuts that are applied to the A(Obj) — 
5'(0bj) plane that were however recovered by the RFF — 
A (Res) diagnostic. This set of galaxies is made of 69 galaxies. 
As a comparison, there are 17 visually classified mergers that 
would have been detected by the usual CAS diagnostic. 

Figure [represents a number of examples taken from the 
69 visually classified galaxies that would have not been re- 
trieved by the A(Obj) — 5'(0bj) diagnostic using its usually 
adopted limits. For each of the selected galaxies, this figure 
presents the COMBO- 17 ID, its environment, its morpho- 
logical type, and its B-band absolute magnitude. Each inset 
also presents an estimate of the contribution from the less lu- 
minous component to the total FLUX_AUTO SEXTRACTOR 
measurement, expressed as a fraction. This is calculated sim- 
ply by dividing the fiux enclosed in the aperture shown by 
the automatic fiux measurement that is performed by SEX- 
TRACTOR. 

Figure [Ml shows that the RFF — A (Res) method pre- 
sented here has detected galaxies that do present inhomo- 
geneities in their light distribution. These inhomogeneities 
typically amount to 10% of the total fiux received from the 
parent galaxy, as estimated by the crude but robust fiux esti- 
mate given by SEXTRACTOR. This confirms the potential 
validity of the approach presented in the current paper to 
separate minor mergers from the whole parent population of 
galaxies. However, Figure [14] also makes manifest the reason 
of the fairly high contamination by non-mergers that has 
been found in ^5.11 It is very difficult, even with the eye, to 
tell apart galaxies undergoing a minor merger episode from 
galaxies that are experiencing other local phenomena such 
as star formation in HII regions. This is particularly true in 
the case of very late minor mergers, in which the less mas- 
sive galaxies have been almost entirely dissolved amidst the 
larger galaxies. However, for less evolved merger systems for 
which the less luminous object has not entirely lost its indi- 
viduality, it is easier to separate mergers and non- mergers. 
Specifically, objects 20213 and 7479 represent two cases of 
objects selected by the RFF — A (Res) that could indeed be 
star formation enhancements. Finally, Figure [TJ] also shows 
that, although the structural merger diagnostic used here 
can indeed select galaxies with inhomogeneous light distri- 
butions, this method alone cannot identify which light clump 
is to be identified as the potential satellite. This is best seen 
in the case of ID=40654. 



6 CONCLUSIONS 

We present a new structural merger diagnostic geared to- 
wards the structural detection of minor mergers which is 
entirely based in the morphological properties of the residual 
images of galaxies after the subtraction of a smooth Sersic 
model. The new indicator makes use of the asymmetry of 
the residuals and of the Residual Flux Fraction of the fit, 
both calculated over the Kron aperture of the galaxies. This 
diagnostic has been objectively proven to be able of pro- 
ducing merger samples of equal or better statistical quality 
than samples obtained using other well established meth- 
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Figure 14. Visually classified minor merger candidates found by 
the RFF — A(Res) diagnostic but missed by the CAS method 
adopting its usual limits. COMBO-17 IDs, environments, mor- 
phological types and B-band absolute magnitudes from the 
STAGES public catalogue are shown for each galaxy. The panels 
also give an indication of the angular extent of the image insets 
and the approximate fractional contribution of the light included 
in the circle to the total SEXTRACTOR automatic flux measure- 
ment. Objects 20213 and 7479 might be localized star formation 
episodes in irregular galaxies. Also, it is not straightforward to 
identify the putative satellite in object 40654, which presents an 
alternative satellite marked in red. 
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ods based on the morphological properties of the original 
images. In particular, objects with symmetric residuals for 
which the Residual Flux Fraction is larger than 0.2 or ob- 
jects with more asymmetric residuals for which the RFF is 
larger than 0.1 are very good candidates to be mergers. We 
have also found that the Gini index of the residuals could 
also produce merger samples of high statistical purity. In 
this case, objects for which the Gini index of the residual 



Figure 14 - continued Visually classifled minor mergers found 
by the RFF — A(Res) diagnostic but missed by the CAS method 
adopting its usual limits. COMBO-17 IDs, environments, mor- 
phological types and B-band absolute magnitudes from the 
STAGES public catalogue are shown for each galaxy. The panels 
also give an indication of the angular extent of the image insets 
and the approximate fractional contribution of the light included 
in the circle to the total SEXTRACTOR automatic flux mea- 
surement. 
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image calculated within the Kron aperture is higher than 
0.5 are also good merger candidates. 

Using the structural parameters of the residuals and the 
Umits provided by the F — score optimization process shown 
here, we have split the whole population of galaxies into two 
different set of galaxies. The first set, sharing the structural 
trends and properties of the mergers included in the train- 
ing set is shown to contain the majority of major and minor 
mergers. The second set has been shown to be almost com- 
pletely free of mergers, as exposed by the very low negative 
contamination rates. However, given the relative dearth of 
mergers among the general galaxy population and the self- 
imposed goal of detecting the more elusive minor mergers, 
it turns out that the RFF — A (Res) diagnostic introduced 
in this paper works best as a negative merger test. In other 
words, it is very effective at selecting non-merging galaxies. 
In common with all the currently-available automatic meth- 
ods, the sample of both major and minor merger candidates 
selected by our test is heavily contaminated by non-mergers, 
and further steps are needed to produce a clean merger sam- 
ple from the first set of galaxies. Nevertheless, the method- 
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ology introduced in this paper can be very useful when ap- 
phed to the huge datasets provided by modern large area 
surveys such as SDSS or UKIDSS. By visually classifying a 
relatively small and manageable number of galaxies, one can 
derive the resulting cuts and best borders which can then 
be applied to the whole sample. This would produce two 
statistically well-defined sets of galaxies, a smaller one con- 
taining the vast majority of both major and minor mergers, 
and a much larger one almost completely devoid of them. To 
identify bona-fide mergers, only the first set would need to 
be processed further to remove the non mergers. This could 
be done either by visual inspection or by using additional 
information such as colours or 3-D spectra. 

This work also suggests that the use of the structural 
parameters of the residual images could indeed be used as 
a tool to study the properties of minor mergers. We argue 
that this is due to the fact that by removing the bulk of 
the host galaxy light we might able to detect much fainter 
merging galaxies and over much longer timescales. This will 
be further studied by using N-body merger simulations in a 
forthcoming paper, which will eventually provide the way to 
link the structural parameters measured over the HST/ACS 
images with the underlying properties of the observed merg- 
ers. 
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